Gauge fields and interferometry in folded graphene 
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Folded graphene flakes are a natural byproduct of the micromechanical exfoliation process. In 
this Letter we show by a combination of analytical and numerical methods that such systems behave 
as intriguing interferometers due to the interplay between an externally applied magnetic field and 
the gauge field induced by the deformations in the region of the fold. 



Gauge fields modulate the dynamics and interactions 
of electrons in a variety of scales from particle cosmology 
to phase transitions in condensed matter pQ. Graphene 
is a one-atom-thick carbon crystal where charge carri- 
ers behave as massless quasiparticles [2 and where the 
impact of lattice deformations on electrons is equivalent 
to effective gauge fields [2H5]. Graphene behaves as a 
highly stretchable membrane [6], whose elastic deforma- 
tions can be induced in a controlled way [7 . In this 
Letter we explore the most natural setup, a graphene 
fold [8] (see also [9]), where the interplay between an ap- 
plied magnetic field and the effective gauge fields leads to 
pronounced interference effects between chiral electronic 
modes. The device allows for the controlled splitting of 
electronic current paths, opening the way for new multi- 
terminal (Mach-Zehnder) interferometers [10] . 

We consider the setup sketched in Fig. [IJi). It consists 
of a graphene nanoribbon with armchair edges (AGNR), 
which has been folded along a line perpendicular to its 
longitudinal axis, here taken to be along the x direction. 
We assume that the two layers are completely decou- 
pled [8| and that they are connected only through the 
region of the fold. Under these assumptions, current in- 
jected through a lead attached to the bottom layer, say, 
can be extracted from a lead attached to the top layer 
after having crossed the region of the fold. 

The radius of the fold R is determined by the bal- 
ance between the bending rigidity of graphene n and 
the van der Waals attraction 7 between the two lay- 
ers [TTJ [12]. We find (see Supplementary Informa- 
tion) R « « 7 A, where n « 1 eV and 7 « 

0.022 eV x A . The deformation associated with the 
fold also induces a uniaxial strain along the x axis, 
Uxx ~ 7/(A + 2/i) = u w 0.1% where A w 3 eV x A" 2 
and \± ~ 9 eV x A are the elastic Lame coefficients 
of graphene. This strain induces an inter- valley gauge 
field [1 [13], A s x tr = 0, Af = f3u/a « 1.6 x HT 3 A" 1 , 
where a is the carbon-carbon distance (a = ao ~ 1.42 A 
at equilibrium), f3 — — <91og(£)/<91og(a) « 2 — 3, t be- 
ing the hopping between ir orbitals in nearest-neighbor 
carbon atoms (t = to ~ 3 eV at equilibrium). The cur- 
vature of the fold hybridizes ir and a bands, leading to 




FIG. 1: (color online), a) A folded graphene ribbon in a mag- 
netic field, b) A topologically equivalent geometry obtained 
by unfolding the ribbon. In this geometry it is clear that the 
magnetic field has opposite signs in the two layers. In the re- 
gion of the fold (shaded area) two co-propagating snake states 
(red and blue lines), which are spatially-separated due to the 
deformation-induced gauge field, are present, c), d) The low- 
energy dispersion relations in the asymptotic leads [c)] and in 
the region of the fold [d)] are plotted as functions of Bloch 
momentum k (in units of 7r/o, a being the appropriate lattice 
constant) . 



another contribution to the gauge field [13] , A^ yh = 0, 
A^ h = (36^/ShvF) x (a/R) 2 « 7.2 x 10" 3 A" 1 , where 
hvp = 3ta/2 w 6.4 eV x A is the Fermi velocity mul- 
tiplied by H, e n7T = V pp7T /3 + V ppa /2 w 3.0 eV, and 
Vpp7r,Vppa are hoppings between p orbitals in nearest 
carbon atoms with different orientations. The small- 
ness of the fold radius implies that the contribution 
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from orbital hybridization is dominant. A gauge field of 
this magnitude varying over the length of the fold gives 
rise to an effective magnetic field with magnetic length 

^eff ~ [(^ R )/ A y yh } 1/2 ~ 5 nm > SO that ^eff ~ 25 T. 

Let us now consider the folded AGNR in a perpen- 
dicular magnetic field B = Bz [see Fig. [IJi)], strong 
enough to quantize the single-particle spectrum into well- 
resolved [15 Landau levels (LLs), i.e. we assume that 
the magnetic length l& = y/he/(eB) is much smaller 
than the width W of the ribbon. Under these condi- 
tions, current in both layers is carried by edge states, 
which are localized on opposite sides of the sample for 
opposite current directions. Upon entering the region of 
the fold while moving along the transport (x) direction, 
the out-of-plane component of the magnetic field first de- 
creases, then becomes zero, and finally changes sign. 

Since we are neglecting inter-layer hopping, our folded 
AGNR is topologically equivalent to an unfolded ribbon 
in the presence of a magnetic field step, as illustrated 
in Fig. [TJd). Let us now imagine to inject current from 
left to right. In this case, current which is first carried 
by edge states near the bottom edge of the left layer 
will have to be carried, after crossing the fold, by edge 
states on the top edge of the right layer. At low energies 
(more precisely, at energies below the first LL plateau), 
only one edge state is involved in carrying this current. 
In the following we will restrict ourselves to this low- 
energy "one-channel" regime. Fig. [TJ^) shows a zoom of 
the dispersion relation E = E(k x ) of an ideal AGNR in 
the presence of a quantizing magnetic field. Only the first 
few low-energy subbands forming the zeroth and the first 
LLs have been plotted. Note that at sufficiently large 
values of \k x \ these subbands become dispersive giving 
rise to the aforementioned quantum-Hall edge states. 

Within (or in the proximity of) the region of the fold, 
a pair of transverse edge modes which live parallel to the 
axis of the fold (the y axis) are induced by the change of 
sign of the magnetic field. These states are depicted by 
vertical blue and red lines in Fig. [TJd). Their existence 
can be understood within a simple semiclassical picture. 
Similarly to what happens at the sample edges, skipping 
orbits form in each separate layer close to the B\ — B 
interface. These pair of skipping orbits living on the op- 
posite sides of the interface merge together in a pair of 
so-called "snake states" [16]. In an ideal, fiat, B\ — B in- 
terface these two snake states are spatially superimposed 
and centered on the interface (which we define it to be 
the line where the B field vanishes). 

The effect of the gauge field in the region of the fold 
is of paramount importance. It shifts the value of the 
momentum, and, as a result, it changes the position of 
the LL guiding center. The effective magnetic field has 
opposite sign for the two valleys, and is maximal in the 
armchair configuration, while it is zero in the zigzag case. 

To understand better the nature of the snake states 
in the proximity of a B\ — B interface in the pres- 



ence of deformations, we introduce an "auxiliary sys- 
tem" consisting of a B\ — B interface infinitely extended 
along the y direction, centered at x = 0, and of fi- 
nite length w to t along the x direction. We assume that 
a finite deformation-induced gauge field is present for 
\x\ < w/2 <C Wtot- By looking at the spectrum E(k y ) 
(Bloch momentum k y is a good quantum number be- 
cause the auxiliary-system interface is infinitely extended 
along the y direction) and at the eigenstates of this aux- 
iliary system we can deduce some quantitative informa- 
tion about the snake states along the y direction in our 
system, which has instead a finite transverse width. 

Performing analytical calculations based on the contin- 
uum model and numerical ones based on a tight-binding 
Hamiltonian, we find the results plotted in Fig.[l]i). Due 
to the deformation-induced gauge field the band struc- 
ture E(k y ) shows non-equivalent K and K' valleys. At 
low energies there is one dispersive edge state per valley 
[red and blue filled circles in Fig. [lji)]. Both edge states 
live close to the interface and are spatially separated by 
a distance S x given by 
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using the parameters described earlier. [Eq. (TJ applies 
only for magnetic fields B such that S x < ttR.] For real- 
istic values of the magnetic field the separation between 
edge states is of the order of the length ttR of the folded 
region, and the carriers near the fold are valley polarized. 
This separation prevents short range scattering from mix- 
ing the valleys, unlike in other proposals to achieve valley 
polarization in graphene. The number of flux quanta in 
the region between the two channels, using the previous 
parameters, is 



(ttR)WB W[nm]B[T] 



1800 



(2) 



where $0 ~ 4.1 x 10 3 T x nm 2 is the magnetic flux quan- 
tum and we are neglecting a factor proportional to the 
angle between the field and the local orientation of the 
layer [we assume that the fold is not perfectly symmetric, 
see Fig. [lji)]. 

The flux in Eq. Q allows for the operation of the de- 
vice in Fig. |TJl) as an interferometer. Incoming electrons 
propagating along the bottom edge state in the left layer, 
say, can either be transmitted to co-propagating edge 
states on the right layer or be reflected back. In both 
cases, the outgoing edge state is localized on the top edge 
of the sample. Electrons can move from one side to the 
other thanks to the existence of the snake states discussed 
above. Due to the gauge field-induced separation S x the 
two snake states enclose a finite magnetic flux, Eq. 
which, in turn, can result in quantum interference. 

We now turn to present our main numerical results 
for the differential conductance G of a folded AGNR in 
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a strong magnetic field and in the presence of deforma- 
tions in the region of the fold. Our calculations, which 
rely on the Landauer formalism, are based on the tight- 
binding description and on the recursive Green's function 
technique [T7] . 

In Fig. [2] we present data for G as a function of the ap- 
plied magnetic field B. The simulated ribbon has a width 
W ~ 50 nm: for this size the quantum Hall regime is 
reached at B > 30 T. The window of magnetic fields dis- 
played on the horizontal axis in this figure can be shrunk 
at will, provided that W is increased (see Supplementary 
Information). The length w of the folded region, in which 
the gauge field is active, is taken to be w ~ 15 nm. We 
describe the total gauge field A = A str + A hyb as a mod- 
ulation of the hopping, although it is mostly due to the 
hybridization of the tt and a bands (see Supplementary 
Information). This description preserves the symmetries 
and main features of the gauge field, and should not 
change significantly the results. The modulation of the 
hoppings parallel to the x axis, St /to = (t — to)Aoj varies 
between 5% [Fig. ^)] and 10% [Fig.^)]. These values 
are higher than those expected in a realistic fold. They 
allow us to illustrate the combined effect of gauge poten- 
tials and magnetic fields in ribbons of widths amenable to 
numerical analysis. A lower gauge field will give similar 
results in wider ribbons, see Eqs. and 

The most important results of this work are the data 
represented by filled circles in Fig. |2j They clearly show 
an oscillatory behavior of the conductance G as the mag- 
netic field is varied. We emphasize that these oscillations 
are only due to the presence of a gauge field in the re- 
gion of the fold. To make this more evident, in Fig. [2] we 
have also reported the conductance of the same B\ — B 
interface in the absence of deformations (filled squares). 
We can clearly see how the conductance oscillations dis- 
appear when the gauge field is switched off. In this case, 
one can see that G is, in the one-channel regime, always 
equal to the maximum conductance 2e 2 //i, apart from 
possible sharp resonances. This is part of a more gen- 
eral result: due to symmetry properties of the zeroth LL, 
a folded graphene nanoribbon in the absence of a gauge 
field shows always perfect transmission when current is 
carried by the zeroth edge state [18] . 

Some comments related to the data in the presence of 
gauge fields are now in order, i) Roughly, we can state 
that the behavior of the conductance as a function of the 
B field is characterized by two quite distinct regimes. 
We can identify a region of type "I" , approximately com- 
prised between 30 and 50 T, and a region of type "II" 
for B > 50 T. (Note that for B < 30 T one exits the 
one-channel regime: for low fields the first LL energy is 
smaller than the chemical potential chosen in Fig. [2]) In 
region I the oscillations in G are rather clear and reg- 
ular, while in region II the conductance increases and 
then remains essentially large and constant for strong 
magnetic fields (modulo sharp Fano-like resonances), ii) 
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FIG. 2: (color online), (a) Numerical results for the two- 
terminal conductance G (in units of 2e 2 /h) of the device in 
Fig.JTJi) as a function of the applied magnetic field (in T). In 
these simulations W ~ 50 nm and w ~ 15 nm. The chemical 
potential is located at w 180 meV above the Dirac level. The 
data labeled by filled (red) circles refer to the case of a folded 
AGNR with a modified hopping integral t/to = 0.95 along 
the x axis in the region of the fold. The data labeled by filled 
(black) squares refer to the case of a folded AGNR in the 
absence of deformations (t/to = 1.0 everywhere), (b): Same 
as in a) but the data labeled by filled (red) circles have now 
been obtained for a value of St /to twice as big, i.e. t/to = 0.90. 



Comparing the data in panel b) with those in panel a), 
we notice that a larger gauge field produces more regu- 
lar and well-defined oscillations, but with smaller ampli- 
tude. Moreover, for larger gauge fields the type-II be- 
havior is reached at larger values of the magnetic field. 
Let us now examine more carefully the nature of the 
eigenstates of the auxiliary system introduced above. In 
Fig.[3]we present numerical results for its eigenfunctions. 
Fig. [3k) shows a low-energy zoom of the dispersion re- 
lation E(k y ) as a function of the Bloch momentum k y 
parallel to the B\ — B interface. We have selected only 
the Brillouin-zone portions related to the snake states we 
are interested in. We immediately observe that the ex- 
istence of two adjacent regions with opposite magnetic 
fields modifies the usual LL structure of an ideal ribbon 
[Fig. by introducing new dispersive branches which 
correspond to snake modes propagating along the inter- 
face. The presence of the gauge field has a crucial role 
in breaking the symmetry between the two valleys, as it 
generates a downward bending in the first LL subband of 
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FIG. 3: (color online), a) A zoom of the dispersion relation 
of the auxiliary system (see main text) in the vicinity of the 
Dirac points, b) 2D color plot of the K- valley squared wave- 
functions \ipk y (x)\ 2 of the auxiliary system (only the value of 
the wave- function on one type of sublattice, say the A type, is 
shown) . In the horizontal axis the spatial position x is shown 
in units of the magnetic length is- In the vertical axis we plot 
the Bloch momentum k y parallel to the interface, c) Same as 
in b), but for the K' valley, d) Same as in the previous panels, 
but for the first-excited subband, which includes the first LL 
and the downward dip induced by the gauge field (only the 
K valley is considered, the K' valley has no dip). 



detailed analysis of the splitting of the LLs near the fold 
is given in the Supplementary Information. 

The results presented here describe two terminal mea- 
surements in a folded graphene layer, and they assume 
that each layer can be contacted separately. The anal- 
ysis can be extended in a straightforward way to multi- 
terminal setups, where more complex correlations can be 
monitored. The combination of flexibility and stiffness 
of graphene allows for many other combinations of gauge 
fields and electronic currents. 
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SUPPLEMENTARY INFORMATION 

This section contains technical details and extra nu- 
merical results relevant to the main text. 



the K valley, say, without modifying the same LL in the 
K' valley. Figs. (3}})-d) illustrate two-dimensional color 
plots of the envelope-function profile ^k y (x) as a func- 
tion of the Bloch momentum k y and of the coordinate 
x along the transport direction, for the three types of 
snake states we are considering. More precisely, Fig. J3}d) 
shows the wave-function belonging to the K' valley. For 
0.55 < kya/ir < 0.75 we have a standard zero-LL edge 
wave-function whose guiding center moves from the right 
side of the auxiliary system toward the center. For larger 
values of k y the wave-function becomes a snake state lo- 
calized close to the interface. In Fig. [3J3) we consider 
instead the K valley, and we plot the wave-function pro- 
file corresponding to the lowest-energy subband. Again, 
for small values of k y we find ordinary edge states with a 
guiding center that moves from the left side of the aux- 
iliary system towards the center. This is then converted 
into a snake state which approaches the interface asymp- 
totically. It is crucial to observe that the "centers" of 
the two snake states corresponding to the two valleys 
are not superimposed but are instead displaced by a fi- 
nite amount. Finally, in Fig. [3]i) we have plotted the 
wave-function ^ (x) of the states belonging to the first 
LL. Since this is an excited subband, the correspond- 
ing wave-functions have a node, and thus the probability 
amplitude \ijjk y (x)\ 2 presents a double-peaked structure 
in space, centered around the B\ — B interface. A more 



Elastic strains and gauge fields 

The elastic energy of a graphene layer near a planar 
configuration can be written as [19] 



J ^ r { ^ l u xx(r) + u yy (r)} 2 + 
where is the deformation tensor defined by [T9] 



^yy &yUy 



(d x u z ) 2 
2 

(dyU z ) 2 



^x^y H - &yU x D x U z 0yU z 



(4) 



Ui being atomic displacement vectors. The parameters 
A « 3 eV x A and [i « 9 eV x A are the Lame 
coefficients of graphene, while k ~ 1 eV is the bending 
rigidity [20]. The van der Waals attraction energy be- 
tween the graphene layers is 



^vdw = 7 / d 2 r 
Jn 



(5) 
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where Q is the contact region, which does not include 
the fold, and 7 « 0.022 eV x A is the van der Waals 
interaction [2T] . 

To a first approximation, the radius of the fold R is 
determined by the balance between the bending rigidity 
and the van der Waals interaction [11]. Approximating 
the fold by half a cylinder of radius R we find 



R, 



7A. 



(6) 



The presence of a substrate and the finite separation be- 
tween the layers, d « 3.3 A, tend to increase this value 
and make the shape asymmetric. The distribution of 
curvatures and strains (see below) will become inhomo- 
geneous and asymmetric as well. 

Now, the bending energy is reduced by an expansion 
of the half cylinder, which leads to in-plane strains. As- 
suming a constant uniaxial strain, u, we find 



E e i(R) 



ttRk ttR(X + 2fi)u 2 



R 2 (l + u) 



(7) 



Minimizing this expression with respect to u and using 
Eq. (pi), we finally find 
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(A + 2fi)R 2 (A + 2n) 



1.0 x 10" 



(8) 



A similar result has been found for carbon nanotubes [22] . 
In the folded region, this strain leads to a constant gauge 
field [3] [131 [23] 



r A!, tr = 



A f = POL « 1.6 x lO- 3 !" 1 
y a 



(9) 



where a is the distance between nearest neighbors in 
the graphene lattice (a = olq « 1.42 A at equilibrium), 
j3 = —d\og(t)/d\og(a) ^2 — 3, and t is the hopping be- 
tween 7r orbitals in nearest neighbor atoms (t = to ~ 3 eV 
at equilibrium). This situation has been considered in 
Ref. |24] In the absence of a magnetic field, the transmis- 
sion can be calculated analytically. 

The bending of the layer induces the hybridization of n 
and a orbitals, and the emergence of a new contribution 
to the total gauge field, not related to strains [14 . Its 
value is 



f Af h = 

^hyb 



3e™ a 2 r2 x 10 -3 A -i 



(10) 



Shv F R 2 



where hv-p = 3ta/2 



6.4 eV x A is the Fermi velocity 

- V^rjTr / 3 H - Vr)1 



multiplied by h and e n7r = V pp7T /3 + V ppa /2 w 3.0 eV, 
V pp7r and V ppcr being the two possible hopping integrals 
between p orbitals in neighboring carbon atoms. Due to 



the smallness of the radius of the fold, the contribution 



from the hybridization in Eq. (10) is roughly one order 
of magnitude larger than the one due to strains, given in 
Eq. ([9]). It is equivalent to a relative change of the value 
of the hopping t of 1%. 
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FIG. 4: Wave-functions in the B\ — B interface in the pres- 
ence of deformations. Panel a) Probability density \ipk y (x)\ 2 
for the snake states residing along the B\ — B interface. The 
horizontal axis is a zoom around the center of the auxiliary 
system where the interface is located. Data labeled by filled 
red circles refer to the K valley, while data labeled by filled 
blue squares refers to the K' valley. The two vertical thin 
lines define the deformed region within which the B — »> —B 
transition takes place. We have chosen an intermediate en- 
ergy value in the one-channel regime, while the magnetic field 
is B = 40 Tesla. Panel b) Probability density \ipk y (x)\ 2 for the 
snake state belonging to the first-excited subband [see Fig. 3d) 
in the main text] . Notice that this single eigenstate is charac- 
terized by a double-peak probability distribution, which could 
generate interference effects similar to the ones occurring in 
an Aharonov-Bohm ring inteferometer. 



Snake states 

Some of the important features that might be hid- 
den in the color plots of Fig. 3 in the main text have 
been highlighted here in Fig. [4j which contains typical 
one-dimensional cuts of ipk y (x) fc> r a fixed value of k y . 
More precisely, Fig. [4^i) reports the squared modulus 
(x)\ 2 of the A-sublattice wave- function correspond- 
ing to a snake state in the lowest-energy subband. The 
spatial region where deformations are present (here taken 
to be 6 nm wide) is inside the area delimited by the two 
vertical lines. The position where the B field vanishes in 
not shown in this plot for simplicity but lies within this 
region. The value of k y chosen to make this plot corre- 
sponds to an energy intermediate between the zero and 
the first LL. As shown in Fig. 3 of the main text, one 
of these two snake states belongs to the K valley, while 
the other one belongs to the K' valley. From Fig. |4^i) 
it is very clear that the centers of the two snake states 
are slightly shifted with respect to each other by roughly 
10 nm. Some finite magnetic flux can thus be enclosed 
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FIG. 5: Conductance through a nanoribbon with zigzag 
edges. Plot of the two-terminal conductance (in units of 
2e 2 /h) of the folded graphene nanoribbon with zigzag edges 
along the transport (x) direction. The parameters used to 
produce this plot are identical to those used in Fig. 2 of the 
main text. For the sake of comparison, the data for armchair 
edges have been reported again. It is evident that in the 
case of zigzag edges no regular and distinct oscillations are 
present, the conductance being almost always equal to unity 
apart from some resonances and minor variations. 



Influence of lattice orientation 

As already shown in Fig. 2 of the main text (data la- 
beled by black filled squares), the conductance oscilla- 
tions discussed above are absent when the gauge field is 
switched off. Our results then are similar to those re- 
ported in Ref. [T8l 

A similar disappearance of the oscillations occurs if 
the orientation of the lattice with respect to the fold axis 
is rotated by 90° (zigzag edges along the transport di- 
rection). The gauge field in this case is A = A(x)x ex 
Q(w/2 — \x\)x. Such a pseudo- vector potential has zero 
curl, i.e. it can be gauged away. We thus expect that 
the oscillations disappear in this case: this is indeed con- 
firmed by our numerical calculations, as illustrated in 
Fig. [5] It is interesting to note that the gauge transfor- 
mation is applicable only to the effect of the gauge field 
on the electronic wavefunctions. The lattice deforma- 
tions which originate the gauge field, obviously, cannot 
be cancelled. 

Numerical results for wider ribbons 



between them, assuming that their spatial displacement 
is not perfectly symmetric with respect to the position 
where B vanishes (this is necessary to avoid a perfect - 
and unrealistic - cancellation between the fluxes coming 
from the B and — B regions). This magnetic flux is at 
the origin of interference effects and conductance oscilla- 
tions shown in Fig. 2 of the main text. If the deformed 
B\ — B interface provides a finite coupling between the 
two snake states belonging to opposite valleys, the en- 
closed residual flux manifests itself as a phase difference 
between the two quantum trajectories involving the K 
and K' interface snake states. We believe that this is the 
mechanism responsible for the conductance oscillations 
of type II shown in Fig. 2 of the main text. 

Another appealing possibility of interference is offered 
by the snake state depicted in Fig. [4)3), which belongs 
to the first subband of the K valley. As we have al- 
ready pointed out, this state presents a node close to the 
interface, and thus the corresponding probability den- 
sity \ifjk y (x)\ 2 possesses a double-peak structure. In this 
case thus there is no need of an inter-valley coupling 
mechanism to observe interference. An incoming elec- 
tron can split at the interface into a two-path configu- 
ration, which can generate Aharonov-Bohm-type inter- 
ference analogous to the one seen in ordinary quantum 
rings. We have checked that the type-I oscillations in Fig. 
2 of the main text occur precisely in the magnetic-field 
range where this excited snake state exists at the chosen 
chemical potential. 
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FIG. 6: Gauge-field-induced conductance oscillations through 
a wider ribbon. Plot of the two-terminal conductance of an 
armchain-folded-graphene nanoribbon with W = 100 nm. 
The region where deformations are present is 30 nm long, 
the chemical potential has been set to ~ 100 meV above the 
Dirac level, and the hopping variation is St /to = 3% [panel 
a)] and 10% [panel b)]. Color and symbol coding are identical 
to those used in Fig. 2 of the main text. 
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As already mentioned in the main text, most of our 
numerical calculations have been carried out in graphene 
ribbons with width W = 50 nm solely for computational 
reasons. The corresponding quantizing magnetic field, 
given by the requirement £b ^ W, turned out to be B > 
30 Tesla. The range of magnetic fields where oscillations 
in the conductance are seen (see Fig. 2 in the main text) 
is largely out of experimental reach. Nevertheless, the 
same analysis can be carried out for wider ribbons. Here 
we have chosen a ribbon with a transverse width of W = 
100 nm for which the quantum Hall regime is reached 
at 10 Tesla, the conductance oscillations occuring in a 
range of magnetic field values (10 < B < 30 Tesla) that 
is currently experimentally accessible. In Fig.|6]we report 
conductance data for this system. The value of St /to is 
3% in panel a) and 10% in panel b), while the length w of 
the folded region is taken to be w — 30 nm. We clearly see 
that Fig. [6] shows a pattern of oscillations which shares 
the same basic physical features with the data reported 
in Fig. 2 of the main text. 
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